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Abstract 

We briefly recall the main physical features of the parton distributions 
in the quantum statistical picture of the nucleon. Some predictions from a 
next-to-leading order QCD analysis are successfully compared to recent un- 
polarized and polarized experimental results. We will discuss the extension 
to the transverse momentum dependence of the parton distributions and its 
relevance for semiinclusive deep inelastic scattering. Finally, we will present 
some new positivity constraints for spin observables and their implications 
for parton distributions. 
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1 A short review on the statistical approach 



Let us first recall some of the basic ingredients for building up the parton 
distribution functions (PDF) in the statistical approach, as oppose to the 
standard polynomial type parametrizations, based on Regge theory at low x 
and counting rules at large x. The fermion distributions are expressed by the 
sum of two terms p], the first one, a quasi Fermi-Dirac function, for a given 
helicity and flavor, and the second one, a flavor and helicity independent 
diffractive contribution equal for light quarks. So we have, at the input 
energy scale Ql = 4GeV 2 , 



AX% q x b Ax b 
exp[(a; — Xq)/x] + 1 exp(x/x) + 1 



Ql) = u ^Iw-t ■ ; + , (!) 



- h( n 2, j^ogVg , Ax b 

xq (x, Qn) = - — 7 1 ; — — . 2 

exp [{x + X^)/x\ + 1 exp(x/x) + l 



Notice the change of sign of the potentials and helicity for the antiquarks. 
The parameter x plays the role of a universal temperature and X^ g are the 
two thermodynamical potentials of the quark q, with helicity h — ±. It 
is important to remark that the diffractive contribution occurs only in the 
unpolarized distributions q(x) = q+(x) and it is absent in the valence 

q v (x) — q(x) — q(x) and in the helicity distributions Ao(x) = q+(x) — q~(x) 
(similarly for antiquarks). The eight free parameters^! in Eqs. ( TT1I21 were 
determined at the input scale from the comparison with a selected set of 
very precise unpolarized and polarized Deep Inelastic Scattering (DIS) data 
[1]. They have the following values 

x = 0.09907, b = 0.40962, b = -0.25347, A = 0.08318, (3) 

X + = 0.46128, Xq u = 0.29766, X~ d = 0.30174, X+ = 0.22775 . (4) 
For the gluons we consider the black-body inspired expression 

xG(x,Ql)= * G f_° , (5) 
exp(x/x) — 1 



2 A = 1.74938 and A = 1.90801 are fixed by the following normalization conditions 
u — u = 2, d — d = 1. 
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a quasi Bose-Einstein function, with be = 0.90, the only free parameter □, 
since Aq = 20.53 is determined by the momentum sum rule. We also assume 
that, at the input energy scale, the polarized gluon distribution vanishes, 
so xAG(x, Qq) = 0. For the strange quark distributions, the simple choice 
made in Ref. [1] was greatly improved in Ref. [2]. More recently, new tests 
against experimental (unpolarized and polarized) data turned out to be very 
satisfactory, in particular in hadronic collisions, as reported in Refs. [3J Hj. 
An interesting point concerns the behavior of the ratio d(x)/u(x), which de- 
pends on the mathematical properties of the ratio of two Fermi-Dirac factors, 
outside the region dominated by the diffractive contribution. So for x > 0.1, 
this ratio is expected to decrease faster for X^ d — x < x < Xq u + x and then 
above, for x > 0.6, it flattens out. 

This change of slope is clearly visible in Fig. 1 (Left), with a very little Q 2 
dependence. Note that our prediction for the large x behavior, differs from 
most of the current literature, namely d(x)/u(x) — > for x — > 1, but we find 
d(x)/u(x) — > 0.16 near the value 1/5, a prediction originally formulated in 
Ref. [5]. This is a very challenging question, since the very high-x region 
remains poorly known. To continue our tests of the unpolarized parton dis- 
tributions, we must come back to the important question of the flavor asym- 
metry of the light antiquarks. Our determination of u(x, Q 2 ) and d(x, Q 2 ) is 
perfectly consistent with the violation of the Gottfried sum rule, for which 
we found the value Iq = 0.2493 for Q 2 = 4GeV 2 . Nevertheless there remains 
an open problem with the x distribution of the ratio d(x)/u(x) for x > 0.2. 
According to the Pauli principle, this ratio is expected to remain above 1 for 
any value of x. However, the E866/NuSea Collaboration [6] has released the 
final results corresponding to the analysis of their full data set of Drell-Yan 
yields from an 800 GeV/c proton beam on hydrogen and deuterium targets 
and they obtain the ratio, for Q 2 = 54GeV 2 , d(x)/u(x) shown in Fig. 1 
(Right). Although the errors are rather large in the high- a; region, the sta- 
tistical approach disagrees with the trend of the data. Clearly by increasing 
the number of free parameters, it is possible to build up a scenario which 
leads to the drop off of this ratio for x > 0.2. For example this was achieved 
in Ref. [7j, as shown by the dashed curve in Fig. 1 (Right). There is no such 
freedom in the statistical approach, since quark and antiquark distributions 

3 In Ref. pQ we were assuming that, for very small x, xG(x, Qq) has the same behavior 
as xq(x, <3o), so we took be = 1 + b. However this choice leads to a too much rapid rise 
of the gluon distribution, compared to its recent determination from HERA data, which 
requires be — 0.90. 
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are strongly related. On the experimental side, there are now new opportuni- 
ties for extending the d(x)/u(x) measurement to larger x up to x = 0.7, with 
the upcoming E906 experiment at the 120 GeV Main Injector at Fermilab 
[8] and a proposed experiment at the new 30-50 GeV proton accelerator at 
J-PARC 0. 

Analogous considerations can be made for the corresponding helicity distri- 
butions, whose most recent determinations are shown in Fig. 2 (Left). By 
using a similar argument as above, the ratio Au(x)/u(x) is predicted to have 
a rather fast increase in the x range [Xq u — x, Xq u + x) and a smoother be- 
haviour above, while Ad(x)/d(x), which is negative, has a fast decrease in the 
x range (X^ — x, X^ d + x) and a smooth one above. This is exactly the trends 
displayed in Fig. 2 (Right) and our predictions are in perfect agreement with 
the accurate high-x data. We note the behavior near x — 1, another typical 
property of the statistical approach, is also at variance with predictions of the 
current literature. The fact that Au(x) is more concentrated in the higher x 
region than Ad(x), accounts for the change of sign of gi(x), which becomes 
positive for x > 0.5, as first observed at Jefferson Lab [T2"] . 
Concerning the light antiquark helicity distributions, the statistical approach 
imposes a strong relationship to the corresponding quark helicity distribu- 
tions. In particular, it predicts Au(x) > and Ad(x) < 0, with almost 
the same magnitude, in contrast with the simplifying assumption Ati(z) = 
Ad(x), often adopted in the literature. According to the COMPASS experi- 
ment at CERN [T3], Au(x) + Ad(x) ~ 0, in agreement with our prediction. 

2 The TMD extension 

We now turn to another important aspect of the statistical PDF and very 
briefly discuss a new version of the extension to the transverse momentum 
dependence (TMD). In Eqs. ([HE]) the multiplicative factors X% q and (X^)- 1 
in the numerators of the non-diffractive parts of g's and if s distributions, im- 
ply a modification of the quantum statistical form, we were led to propose in 
order to agree with experimental data. The presence of these multiplicative 
factors was justified in our earlier attempt to generate the TMD [13], but it 
was not properly done and a considerable improvement was achieved recently 
[15] . We have introduced some thermodynamical potentials Yq q , associated 
to the quark transverse momentum kj<, and related to X(} q by the simple 
relation ln(l + exp[F ^]) = kXjj . We were led to choose k = 3.05 and this 
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method involves another parameter p 2 , which plays the role of the tempera- 
ture for the transverse degrees of freedom and whose value was determined 
by the transverse energy sum rule. We have calculated the pr dependence of 
semiinclusive DIS cross sections and double longitudinal spin asymmetries, 
taking into account the effects of the Melosh-Wigner rotation, for tt^ pro- 
duction by using this set of TMD statistical parton distributions and another 
set coming from the relativistic covariant approach [16J. Both sets do not 
satisfy the usual factorization assumption of the dependence in x and fcr and 
they lead to different results, which can be compared to recent experimental 
data from CLAS at JLab, as shown on Figs. 3. 

3 Positivity bounds 

Spin observables for any particle reaction, contain some unique information 
which allow a deeper understanding of the nature of the underlying dynam- 
ics and this is very usefull to check the validity of theoretical assumptions. 
We emphasize the relevance of positivity in spin physics, which puts non- 
trivial model independent constraints on spin observables. If one, two or 
several observables are measured, the constraints can help to decide which 
new observable will provide the best improvement of knowledge. Different 
methods can be used to establish these constraints and they have been pre- 
sented together with many interesting cases in a review article [IB] . For lack 
of space, here we will only briefly discuss some new results obtained very 
recently [HI [20]. 

Let us consider the inclusive reaction of the type A (spin 1/2) +B (spin 1/2) —> 
C + X, where both initial spin 1/2 particles can be in any possible directions 
and no polarization is observed in the final state. The spin- dependent cor- 
responding cross section a(P a ,Pb) = Tr(Mp), can be defined through the 
4x4 cross section matrix M and the spin density matrix p, where P a , Pb are 
the spin unit vectors of A and B, p = p a ® pb is the spin density matrix with 
Pa — {I2 + Pa ' <?a)/2, and similar for p&. Here I2 is the 2x2 unit matrix, and 
= (°s> a yi a z) stands for the 2x2 Pauli matrices. M can be parametrized in 
terms of 8 parity-conserving asymmetries and 8 parity-violating asymmetries. 
The crucial point is that M is a Hermitian and positive matrix and this allows 
to derive some positivity conditions. Since one of the necessary conditions 
for a Hermitian matrix to be positive definite is that all the diagonal matrix 
elements has to be positive Ma > 0, we thus derive 1 ± A NN > \A aN ± Aw\, 
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valid in full generality, for both parity-conserving and parity-violating pro- 
cesses, where A N denotes the single transverse spin asymmetry and A NN the 
double transverse spin asymmetry. In the case p^ + p^ — > C + X where the 
initial particles are identical, we have A a N(y) = — y). Using this rela- 
tion , one obtains, 1 ± Ann{v) > |^4.jv(2/) -F An(—v)\- This is an interesting 
result which, can be used, in principle, with available data on An for tt , 

7r°, T] production, to put some non trivial contraints on A NN (y). 
Let us now study the implications of the above relation for the parity- 
violating process p^ + p^ — > W ± + X. Since Ann ~ 0, to a very good 
approximation, it reduces to 1/2 > \A^{y = 0)|, to be compared with the 
usual trivial bound 1 > \A N (y = 0)|. 

The TMD quark distribution in a transversely polarized hadron can be 

expanded as f q/h t(x, k ± , S) = f q /h(x,k±) + ±A N f q/hr (x, k±) S ■ (p x kj, 

where p and are the unit vectors of p and k^, respectively f q /h(x, k±) is 
the spin-averaged TMD distribution, and A N f q / h ^(x, k±) is the Sivers func- 
tion. There is a trivial positivity bound for the Sivers functions which reads 
\A N f q /ht(x, k±)\ < 2f q /h(x,k±). Since A N is directly expressed in terms of 
A JV f q /ht(x, k±), this trivial bound can be improved as shown in Ref. [20] . 
In the helicity basis it is easy to obtain the explicit form of M and now from 
Mjj > 0, we have 1±All(v) > \A a L(y)±AbL,(y)\, where Al denotes the single 
helicity asymmetry and All the double double asymmetry. It is important to 
note that for identical initial particles scattering, one has A aL (y) = A bL (—y), 
so one gets 1 ±All(v) > ± A L (—y)\. These bounds should be tested 

in RHIC experiments for W ± or Z° production in longitudinal pp collisions, 
pp — > W ± /Z° + X. In perturbative QCD formalism, at leading-order and 
restricting to only up and down quarks, one has simple expressions for the 
single and double helicity asymmetries, involving only quark helicity distri- 
butions. The statistical PDF satisfy the positivity bound. Finally at y = 0, 
since All(0) is expected to be very small, the bound implies ^4^(0) < 1/2, 
a remarquable simple result which must be satisfied by future experimental 
data. 
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Figure 1: Left : The ratio d(x)/u(x) as function of x for Q 2 = 4GeV 2 (solid 
line) and Q 2 = lOOGeV 2 (dashed-dotted line). Right : Comparison of the 
data on d/u(x,Q 2 ) from E866/NuSea at Q 2 = 54GeV 2 0, with the predic- 
tion of the statistical model (solid curve) and the set 1 of the parametrization 
proposed in Ref. [7] (dashed curve). 
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Figure 2: Left : Quark and antiquark helicity distributions as a function of 
x for Q 2 = 3GeV 2 . Data from COMPASS [10]. The curves are predictions 
from the statistical approach. Right : Ratios (Am + Am) /{u + u) and (Ad + 
Ad)/(d + d) as a function of x. Data from Hermes for Q 2 = 2.5GeV 2 [H] 
and a JLab Hall A experiment p2]. The curves are predictions from the 
statistical approach 
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Figure 3: Left : The double longitudinal-spin asymmetry Ai for 7r + produc- 
tion on a proton target, versus the tt + momentum p?, compared to the JLab 
data Ref. [T7]. The solid lines are the results from the TMD statistical dis- 
tributions [J5] and the dashed lines correspond to the relativistic covariant 
distributions [TB] Right : Same for 7r~. 
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